Electron Transport in Double Quantum Dot governed by Nuclear Magnetic Field 
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We investigate theoretically electron transfer in a doble dot in a situation where it is governed by 
nuclear magnetic field: This has been recently achieved in experiment We show how to partially 
compensate the nuclear magnetic field to restore Spin Blockade. 
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Much modern research is devoted to practical realiza- 
tion of quantum manipulation and computation (QMC). 
Although QMC can be performed with convenient mag- 
netic resonance techniques 0, this necessary involves 
macroscopically many identical spins. The challenge is 
to do QMC with individual spin states, e.g. those of lo- 
calized electrons. Remarkable experimental progress has 
been recently achieved in preparation, manipulation and 
measurement of individual spin quantum states in quan- 
tum dots. An important issue to resolve in this case is 
the spin measurement: to do this, one has to convert spin 
into charge and/or electric current [^Q. Such conversion 
has been realized in a single quantum dot 0, 0| . Other 
experiments were focused on the transport through two 
coupled quantum dots 0, Hi • Although such double dot 
is a more complicated system with many additional pro- 
cesses influencing spin and charge transfer, the advan- 
tage is an immediate access to spin-charge conversion. 
In a double dot, two electron spins can be entangled 
in the course of quantum manipulation forming either 
symmetric spin singlet or antisymmetric triplet states. 
This strognly affects electron transport giving rise to Spin 
Blockade of electron tunneling The quantum dots are 
commonly fabricated in GaAs-based semiconductor het- 
erostructures. The specifics of GaAs is a strong hyper- 
fine interaction between electron and nuclear spins [10|. 
Therefore, the spin of an electron localized in a a quan- 
tum dot can be strongly affected by the effective spin 
magnetic field B'^ arising from random configuration of 
many nuclear spins situated in the dot. This field helps 
the transitions between the components of spin doublet 
[m as well as between singlet and triplet states It 
has been experimentally proved that the nuclear field not 
only lifts Spin Blockade in a double dot but gives rise 
to time-dependent spin-driven oscillations of the current 
|l3j . Similar effects have been also observed in Quantum 
Hall constriction The origin of the oscillations is 

the modulation of the current by nuclear field and feed- 
back of electron spin on nuclear spins that results in their 
nutations Since the nuclear field is random and 

hardly controllable, its influence on the electron spin sig- 
niflcantly complicates QMC. This has motivated inten- 
sive research aimed to measure and to predict the effect 
of the nuclei on spins in quantum dots l3j la| ■ 

The study presented here has been stimulated by very 



recent experiment in this direction 0. The advantage 
achieved in this experiment is the better control of elec- 
tron levels in the dots and, most importantly, the pos- 
sibility to control and tune the tunnel coupling between 
dots in a wide range. Most interesting results were ob- 
tained near the boundary of Coulomb diamond where 
the states that differ by number of electrons in the dot 
are aligned in energy. The authors of P| were able to 
demonstrate the order-of-value change of the current by 
anomalously small external magnetic field B^xt — 5mT 
that matches the nuclear field. By tuning the tunnel cou- 
pling, they observe this effect in several different regimes. 

We provide an adequate theoretical framework for this 
experimental situation. We derive and solve density ma- 
trix equation valid in the regimes of interest. We con- 
centrate on the fact that nuclear field randomly changes 
at time scale bigger than that of electron dynamics but 
smaller than the measurement time and therefore pro- 
vides "frozen disorder" for electron spins. We achieve 
agreement with experiment in rather fine details. Very 
imortant result of our analysis is that for any given con- 
figuration of nuclear magnetic field there is always a value 
of external magnetic field at which there is no current — 
stopping point. This encourages us to speculate that the 
fast current measurement in this setup can be used to 
measure and, via external feedback, partly compensate 
the nuclear magnetic field. The setup would be stabilized 
in the stopping point where Spin Blockade is resored and 
QMC is possibly enabled. 

The charge configuration of the double dot is given 
by the number of electrons in the left and right dot 
(A^LjTVfl). The gate and bias voltages are tuned in 
the experiment to provide the following transport cycle: 
(1,1) (0,2)^r«(o, l)^rz,^j^]^-)^ rp^Q i^g^ processes 

are incoherent tunnel transitions with electron transfer 
from the left and to the right lead, the tunneling rates 
being and Tn respectively. The first transition may 
be coherent and is due to tunnel coupling t between elec- 
tron states in different dots. If spin is taken into account, 
there are four possible quantum states in (1, 1) config- 
uration: a singlet S'(l,l) and there components of the 
triplet ri(l,l). As to (0,2) configuration, only a singlet 
state 5*2(0,2) participates in the transport: The triplet 
states of (0, 2) are much higher in energy owing to strong 
exchange interaction between two electrons in the same 
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dot. The essence of Spin Blockade is the spin selection 
rule for (1, 1) (0, 2). Provided spin is conserved, there 
is no matrix element connecting any triplet state r(l,l) 
and Sg{0,2). Therefore the transition does not take place, 
the system gets stuck in one of the triplet states and the 
current is blocked. 

The part of the Hamiltonian for (1, 1) and (0, 2) con- 
figurations that conserves spin is presented in the triplet- 
singlet states basis {Ti,S and Sg) as 

Ho = E (^\S){S\+J2\Ti)m\^ + iE + A)\Sg){Sg\ + 

+ti\s){Sg\ + \Sg){sm 

Here A is the detunning of (1, 1) and (0, 2) states linear 
in the gate and bias voltages. The experiments were con- 
centated at the edge of the Coulomb blockade diamond 
where |A| <C eV,Ec- The tunnel coupling between the 
dots mixes two singlets at A ~ i but does not alter triplet 
states. (Fig. ^ 

The leakage current in spin blockade regime can only 
arise from the spin-dependent interactions that mix sin- 
glet and triplet states. Theoretical ly, such interactions 
can be caused by many mechanisms |l6l | . Experimenally, 
the most relevant one appears to be hyperfine interac- 
tion with nuclear spins. Since there are many nuclear 
spins interacting with an electron state in each dot, their 
net effect can be presented in terms of classical vari- 
ables: effective fields B^^.(we measure fields in energy 
units) In the absence of net nuclear polarization, these 
fields are random de pen ding on a concrete configura- 
tion of nuclear spins [llj. Owing to central limit the- 
orem, the distribution of both fields is Gaussian with 
< >= B% = El/N^s, Er, « 0.135meV for GaAs 
being the energy splitting induced by fully polarized nu- 
clei, Ncs being the effective number of nuclei in the dot, 
Ncs ~ 10^~^ for typical dots. It is important for our 
approach that nuclear fields change at time scale of nu- 
clear spin relaxation (~ 1 s), that is much bigger than 
any time scale associated with electron transport. This 
is why they can be regarded as stationary random fields. 
The electron spins inside the dots feel therefore effective 
stationary fields described by 

^^spin = ■ Sl + ■ Sij + Bcxt{Sl + S^) (2) 

Sl,r being the operators of the electron spin in each dot 
and the external magnetic field is || z. We rewrite this in 
triplet-singlet representation as 

i 

2^ \Tn){T±i \ + -j= \S){T±i\ + h.c.j(i) 

where Ba s = (Sf T B^)/2 and 4 = -1,0,1 is the 
projection of the spin of \Ti) state on z-axis. We see 



that the sum of effective fields mixes and splits triplet 
components only. The difference of the fields mixes the 
spin singlet 5'(1,1) and triplet T(l,l) states, this being 
the source of leakage current. 

The energy levels of the resulting Hamiltonian Hgt — 
Ho + ^fspin are determined now not only by the tunneling 
t and misalignment of the levels A but also by the fields, 
the corresponding energy scales can be comparable. The 
mixing of the singlet and triplet in the eigenstates of the 
Hamiltonian can be significant as well. Already from 
analisys of this simple Hamiltonian we can conclude that 
the current is absent if either Bg \\ Ba or Bg -L Ba, 
since here Bg consists of the external and sum of nuclear 
magnetic fields. To see this explicitly from let us 
choose z-axis in the direction of Bg. If Bg \\ Ba, two 
triplet states \T±i) are not mixed with the singlet. If 
Bg _L Ba, it is one state |To) that is not mixed. In both 
cases the system stucks in one of the non-mixed triplet 
states resulting in no current. Importantly, the stopping 
point {Bg, Ba) = can be achieved at any configuration 
of nuclear fields by adjusting i?oxt- 

To evaluate the current in general situation, we pro- 
ceed with formulation of a suitable density matrix ap- 
proach first elaborated for double quantum dot in |l7| . 
Current for the trasport cycle given is proportional to 
the probability to find a system in the state Sg, I = 
erRpSgSg- Although the transport involves 7 states, the 
probabilities of (1,0) doublets are readily expressed via 
other probabilities. So the density matrix to work with 
is spanned by five singlet-triplet states discussed. Using 
the equations of motion, we derive the equations for the 
stationary density matrix (dp/dt = 0). Five diagonal 
equations read 

irflps,s,-*(7^d[^st,/5]|r,) = 

^r^ps^g^ - FinPSS - i(<5'|[^st,p]|S') = 

-r^psgs, + Tinpss - i{Sg\[H,t,P\\Sg) - (4) 

where, motivated by experiment, we also include inelastic 
transitions between S and Sg with the rate Tin , Fm ^ if 
A < 0. The commutator terms include non-diagonal ele- 
ments of density matrix, so we also need 20 non-diagonal 
equations, 

-^(r-''+F>jk-*(j|[i?st,p]|fc)=0. (5) 

Here j,k — Ti,S,Sg number the five states basis,j ^ 
k. The "rates" Fj are zero for triplet states and are 
^Ri^in) for Sg (S*). To close the set of equations we use 
the normalization condition for density matrix, ^ pjj + 
PS,sJI + Tr/Tl)^1,3^T„S. 

The solution gives the current for a given realization 
of nuclear fields. Normally, one expects self-averaging 
over different realizations at time scale of a single mea- 
surement. Since nuclear relaxation times are large, this 
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point deserves some discussion. In fact, raw data aqui- 
sition time in experiment was 0.1 s per point, which 
is probably less than the relaxation time. However, the 
raw data are noisy (see Fig. Q due to both instrumen- 
tal noise and random changes of nuclear fields. An accu- 
rate measurement requires, say, 50 data points, this takes 
time much bigger than the relaxation time. This leads 
us to two conclusions: (i) smooth part of experimental 
data corresponds to the current averaged over realiza- 
tions, (ii) a realistic (factor of 30) improvement of the 
measurement speed and accuracy will allow to measure 
current for a given realization. So that, to compare our 
theory with experimental results, we average the current 
obtained from the solution of equation set 0),© over 
Gaussian distribution of fields. 

Both solving and averaging can be easily done numer- 
ically. To present the physics behind, we give analytical 
results in two limiting cases. The first, natural limit cor- 
responds to small nuclear fieds, Bn <^ max(i, i?cxt)- lu 
this case, the system is preferentially in one of the triplet 
states whose energies are 0, ±Bs. It is convenient in this 
case to choose spin quantization axis along Bs and work 
with parallel and perpendicular components of Ba, Ba'^ 
with respect to this axis. The current reads 



Tfl e/I = 



F{Bs) = t^ + Bl[2 + {Bl + ^^)/f), (6) 

where the first term is due to transitions from \To) and 
the second due to transitions from iTii). As expected, 
the current stops if either b\ = or B^ = 0. The 
average current in this limit 
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We plot the results in this limit for average curent as well 
as for two arbirtaty realizatons of the field (Fig. O . The 
stopping points at -Bcxt — Bj^ are visible for realizations, 
while no features in average current are seen in this range 
The alernative limit of big fields is achived provided 
Bcxt: Bn ^ max(t, t^/A). In this case, the system sticks 
in one of the four states (1,1) with energies ±(5^^ ± Bn). 
The current is determined by decay from these states and 
reads 



I/e = r, 



' 8(A2 + S2) 



(8) 



where n^^^j are unit vectors in the direction of B^ ji. At 
Bcxt = I /e — Tjit^ /B^ and it drops significantly at 
B^^t - max(BN, A). 

Eq. ISI seemingly contradicts to our general statement 
giving non-zero current at Bs J- Ba- A fine point is that 
Eq. ISjis not valid in close vicinity stopping point where 



two of the four states are degenerate. One has to take 
into account that this degerenacy is lifted by coupling to 
Sg. As a result, the current develops a narrow Lorenzian- 
shaped dip in the vicinity of the stopping point. 
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5B — Bext — -Sgl^j, Bw ~ — 



l + i6B/B^y' 
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Since the dip is narrow (see Fig. it is washed away 
upon averaging. 

The average current 
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const iV^w A,B<^Bi, 



/9) 

We encounter a similar situation under conditions 
where the four states are emptied by inelastic tunnel- 
ing. Tin ^ rfl;(i/A)^ while their splitting is determined 
by magnetic field, B^ 3> Ast- It is experimentally con- 
firmed that this always takes place at sufficiently big neg- 
ative A. The current is again an inverse of the sum of 
inverse partial rates and reads 



I/e = Tin {riL X ur) 
The average current as a function of B^xt becomes 



(/)/e = T^nS{Bext/BN); 



(10) 



(11) 



where 



S{x) = 4/a;2 - + V2^erfi(a;/%/2)(6/a;'^ 
-2/x^)exp{-x^/2) - 37rerfi2(a;/\/2)exp(^a;2)/a;^ (12) 

In is interesting to note a special form of this function: 
The graph of S gives a peak with flat top, S"'(0) = 0. 
This funbction provides an excellent fit to experimen- 
tal data (Fig. 0J, those are impossible to fit with more 
conventional peak functions. Such flat peaks are thus 
specific for the model in use and provide strong support 
of its experimental validity. 

In conclusion, we have presented the theoretical frame- 
work for the electron transport via a double quantum dot 
influenced and governed by nuclear magnetic field. Our 
approach is based on density matrix equations and we 
achieve good agreement with experiment Q assuming 
averaging over realizations of nuclear fields. An impor- 
tant feature which is yet to be observed in the course of 
faster and more accurate measurement is the presence of 
stopping points for any given realization of nuclear fields. 

If one interprets the effect of nuclear magnetic fields 
in terms of spin coherence time, the results of are 
discouraging if not forbidding for QMC in GaAs quan- 
tum dot systems. The coherence time estimated is just 
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FIG. 1: Energies of the states with two electrons. The tunnel 
coupling between the dots mixes singlet states S, Sg and does 
not influence triplet states (split by magnetic field). 
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FIG. 2: Average current (solid) and that for two random 
realizations of nuclear magnetic field in the limit of validity of 
^{t/BN = 50, A/t — 1). Note stopping points at _Bext — Bm 
seen for the realizations. 



too short, ~ 10~^ s. We speculate that the presence of 
stopping points can remedy the situation. Faster cur- 
rent measurement would allow to characterize and, with 
the aid of external feedback, partially compensate the 
nuclear fields by stabilizing the system in the stopping 
point. 
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our attention to the topic, many useful! discussions and 
communicating their results prior to publication. We ac- 
knowledge the financial support by FOM. 
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